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Abstract. 

In this letter we show that the force-free Duffing- van der Pol oscillator is completely 
integrable for a specific parametric choice. We derive a general solution for this 
parametric choice. Further, we describe a procedure to construct the transformation 
which removes the time dependent part from the first integral and provide the general 
solution by quadrature. The procedure is shown to have a wider applicability through 
additional examples. We also show that through our method one can deduce linearizing 
transformations in a simple and straightforward way and illustrate it with a specific 
example. 



1. Introduction 

One of the well-studied but still challenging equations in nonlinear dynamics is the 
Duffing- van der Pol oscillator equation [1. Its autonomous version (force- free) is 

x + (a + (3x 2 )x — jx + x 3 = 0, (1) 

where over dot denotes differentiation with respect to time and a, (3 and 7 are arbitrary 
parameters. Equation (JTJ) arises in a model describing the propagation of voltage pulses 
along a neuronal axon and has received a lot of attention recently by many authors. 
A vast amount of literature exists on this equation, for details see for example PJ E] 
and references therein. In most of the works non-integrability properties have been 
considered primarily since the forced version of equation (jlj exhibits a rich variety of 
bifurcations and chaotic phenomena. When (3 = 0, equation becomes the (force- free) 
Duffing oscillator whose integrability and non-integrability properties are well-known pQ, 
while when the cubic term is absent it is the standard van der Pol equation. 

As far as the integrability properties of equation (0) are concerned not much 
progress has been made mainly due to the fact that it does not pass the Painleve 
test as it admits a movable algebraic branch point and a local Laurent expansion in the 
form 3j 
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where r = (t — 1 ) and to an d 03 are arbitrary constants. However, it is known that the 
system (0) admits nontrivial symmetries and a first integral for a specific parametric 
choice [4J 

« = l = ~jr (3) 
The associated first integral reads 

ei* + -a; + ^x 3 ^j = I. (4) 

Rewriting (jH) we get 

+ ^x 3 = Ie~^\ (5) 

which is nothing but a special case of the Abel's equation of the first kind j^j. However, 
it has not been explicitly integrated directly due to the explicit time dependent part 
in the first integral and so the complete integrability of for the choice Q remains 
unclear till date. In this work, we integrate the equation explicitly for the parametric 
choice (JHJ) and provide a general solution for the above specific parametric choice of the 
Duffing-van der Pol nonlinear oscillator, thereby establishing its complete integrability. 
Further, we show that our procedure has a wider applicability, including linearization, 
through additional examples. 



2. Integration of Duffing- van der Pol oscillator for the specific parametric 
choice flU) 

Rewriting {I]) with the specific parametric choice Q we get 

x+(^ + (3x 2 ^jx+^x + x 3 = 0. (6) 

Now introducing a transformation 

w = — xe? 1 , z = e~^*, (7) 

where w and z are new dependent and independent variables respectively, equation (jUJ) 
can be transformed to 

w" - ^-w 2 w' = 0, (8) 

where prime denotes differentiation with respect to z. Equation (JBJ) can be integrated 
trivially to yield 

w' - ^w 3 = /, (9) 
o 

where I is the integration constant. Equivalently, the transformation (J2J) reduces (J3J) to 
this form. Solving 0, we obtain 
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and z is the second integration constant (see figure 1 for a plot of the 
solution x(t)). Replacing the second term inside the square bracket in (jlUJl in terms of 



logarithmic function we get 



i^/~3 log 
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where w = — . This can be further simplified to 
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Rewriting w and z in terms of old variables one can get the explicit solution for the 
equation ©• 




Figure 1. The solution plot of the Duffing- van der Pol oscillator equation JfJJ in the 
form when (3 = 3. 

One may note that a simple solution for equation (JHJ) can also be built from the 
first integral (JSJ). Restricting / = in equation (0) yields 



\/3 . , i 
w = — (2o - z) 2. 



(13) 



Rewriting (fTB^) in terms of the old coordinates we get a particular solution for (0) in the 
form 

V3 



x 



to > 1. 
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3. Method of constructing transformation 



One can construct the transformation (J2J) systematically from the first integral (jH). For 
example, let us split the functional form of the first integral / into two terms such that 
one involves all the variables (t, x, x) while the other excludes x, that is, 

I = F 1 (t,x,x) + F 2 (t,x). (15) 
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Now let us split the function F\ further in terms of two functions such that Fx itself 
is a function of the product of the two functions, say, a perfect different iable function 



d G\(t,x) and another function G 2 (t, x, x), that is, 



at 



1 = F \G^)Jt G ^ x) ) +F ^ x)) - (16) 

We note that while rewriting equation (jl5j) in the form (|16[). we demand the function 
F 2 (t, x) in (fTo"]) automatically to be a function of Gi(£, x). Now identifying the function 
G\ as the new dependent variable and the integral of G 2 over time as the new 
independent variable, that is, 

w = Gi(t, x), 

z = G 2 (t' , x, x)dt' , (17) 



one indeed obtains an explicit transformation to remove the time dependent part in 
the first integral (@J). We note here that the integration on the right hand side of (fTTj) 
leading to z can be performed provided the function G 2 is an exact derivative of t, that 
is, G 2 = ^z(t,x) = xz x + zt, so that z turns out to be a function t and x alone. In 
terms of the new variables, equation (j!6J) can be modified to the form 

-) + fi W . 

\=I-F 2 (w). (18) 




dw\ 
dz) 

Now rewriting equation (fTSjl one obtains a separable equation 

£-/(•), (19) 
dz 

which can be integrated by quadrature. 

For the special case of the Duffing- van der Pol oscillator equation (JB|L by following 
the procedure given above, we can deduce the transformation flJJ) systematically from 
the first integral (@J). For this purpose, we rewrite equation (0) as 




2 



e^ + — x 3 e§\ (20) 
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where J = f 7. Equation (|2*0)l can be further split into 

— e/ 3 — -iC H — 
2 V / 6 

Comparing equation (j2T|) with (fTH|) we get 

Gi = -xe^\ G 2 = |e"^, (22) 

which in turn provides the transformation coordinates, w and z, through the relation 
(fTTj) of the form 
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4. Applications 

We observe that the above procedure can be extended to solve a class of equations. In the 
following we briefly describe the applicability of this method to Duffing oscillator and the 
equation describing the motion of a gaseous general-relativistic fluid sphere. A detailed 
connection between the integrating factors, integrals of motion and the transformation 
coordinates will be published elsewhere 



4-1- Duffing oscillator 

Let us consider the Duffing oscillator [l j 

x + ax + fix + x 3 = jcosut, (23) 

where a, fi, 7 and to are arbitrary parameters. In the absence of the external forcing 
(7 = 0), equation ()23|) has been shown to be integrable for a specific parametric choice 
[Zj, namely, 2a 2 = 9fi. The first integral for this parametric choice has been constructed 
in the form 

r • 2 • 2 2 4 t 

i n Ax axx ax r 

I = e 3 at 1 1 1 . (24) 

[ 2 3 18 4 J v ' 

Equation (|24|) is a complicated nonlinear first order ordinary differential equation (ODE) 

with explicit time dependent coefficients and so cannot be integrated directly. A usual 

way to overcome this problem is to transform equation ()24j) to an autonomous equation 

and integrate it. Using our procedure one can construct the required transformations 

systematically. For example, rewriting equation ()24j) in the form Q15j) we get 

1 / ax\ 2 4 , x 4 4„, 
J=-f± + — ) es at + — e* at . (25) 

Now splitting the first term in equation (J23j) further in the form of (|16|). 

2 / 
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and identifying the dependent and independent variables from ()26|) using the relations 
(|17|) . we obtain the transformation 

X at 3 at_ 

w = —^xe 3 z = e 3 . (27) 

y/2 a K J 

One can easily check that equation (J24j) can be transformed into the autonomous form 

with the help of the transformation (}2Tj) and the latter can be integrated in terms of 

Jacobian elliptic function [7J. We note that the transformation (|2*7)l exactly coincides 

with the one known in the literature which has been constructed in an ad hoc way. 

4-2. Static gaseous general relativistic fluid sphere 

Recently, Duarte etal [S] have constructed the first integral for an equation describing 
the static gaseous general relativistic fluid sphere P, 
t 2 x 2 + x 2 - 1 
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through the so called Prelle-Singer procedure, in the form 
2txx + x 2 + t 2 x 2 — 1 



(29) 



t 2 x 2 

We note that unlike the earlier two examples, namely Duffing and Duffing- van der Pol 
oscillators, we have no exponential function in the first integral. Also the first integral 
is a rational function. However, in the following we show that one can integrate it and 
obtain a general solution for this problem. By following our procedure, equation (|2~H|) 
can be recast in the form 

( tx ) t 2 x 2 ^ ^ 

One may note that the first term in ()30|) can be written in terms of a perfect differential 
form, that is, (j^[log(£x)]) 2 so that G\ = \og(tx) and the function G 2 turns out to be a 
constant, G2 — 1, in the present example. As a consequence we obtain a transformation 

w = \og(tx), z = t, (31) 

which can be utilized to rewrite equation (f2H|) in the autonomous form 

w' 2 - exp(-2w) = I. (32) 

Equation ()32|) can be integrated under two different choices jTOj, namely, > and 
< 0. The respective solutions are 

f\ 1 1 f Vl + e- 2 ™-^ \ 

W Z - Z0 = -WJ lOg { vi + e^ + Vl h />0 ' 
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{ii) z — z = — — — _ arctan y — — — _ — j , / < 0, (33) 
where zq is the second integration constant. Rewriting equation (J3~3|) we get 

{/ 1 _ e -2y"l(z-z )\ 
l0g \2jle-^-*°) J' />0 ' (34) 
lcg( cos ^gp°) ), J<0. 

Utilizing (|3*T|) in one can write down the solution for the static gaseous general 
relativistic fluid sphere of the form 

-^sinhV7(i-t ), I > 0, 
a; = <j v /* (35) 

-cos V-I(t - t ) , / < 0. 



To our knowledge this solution is new to this problem. 



5. Linearization 



Interestingly, our method not only helps to integrate the first integral and gives the 
complete solution for a class of second order nonlinear ODEs but it also helps to deduce 
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the transformations which can be used to linearize the given nonlinear ODEs in a very 
systematic way. For example, let us consider the modified Emden equation, 

x + axx + /3x 3 = 0, (36) 

where a and (3 are arbitrary parameters. Mahomed and Leach JT] have shown that 
equation (|36|) is linearizable through point transformations for a particular parametric 
choice, namely, a 2 = 9/3. The linearizing transformation and the first integral are known 
to be 

t at 2 a 1 

w = — , z=-t--, (37) 

x o 3 x 

and 

I = -t + 7^—r- (38) 

/vail I ry* 

Substituting (J37j) into (|3l)|) one can transform the latter into the free particle equation, 
that is, 

d 2 w , , 

d? = a (39) 

The authors have derived the linearizing transformation through Lie symmetry analysis 
However, in the following we show that the linearizing transformation (J37)l can 
also be derived from the first integral itself through our method. 
Rewriting the first integral (|3*H|) in the form 

d ( t at 2 ' 



x 2 



Q; rr* 2 I rr* 
^ I «AV 



dt \x 6 
and identifying P0|) with (fTH|) we get 



(40) 



t at 2 ^ a x 



Gi = -, G 2 = - + -, F 2 = 0. (41) 

x o 3 

One can see that, unlike the earlier examples discussed so far, in the present example 
G*2 depends both on x and x. However, as we noted earlier the new variable z can be 
obtained once the function G 2 is a perfect derivative of t, that is, 
^ a x d (a i 1\ 

so that (|T7|) gives 

t at 2 1 

x 6 ' 3 x' 

which is nothing but the linearizing transformation. One may note that in this case 
while rewriting the first integral / (equation (jHEI)) in the form (fl3j) . the function F 2 
disappears, that is F 2 = 0, and as a consequence we arrive at (vide equation (fTHjl ) 

^ = / (43) 

which in turn gives (J39)) by differentiation or leads to the solution by an integration. On 
the other hand vanishing of the function F 2 in this analysis is precisely the condition 
for the system to be transformed into the free particle equation jB]. 
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6. Conclusion 

One of the common problems in solving second order nonlinear ODEs is how to solve 
the time dependent first integral and obtain a general solution associated with the given 
equation. In fact most of the methods available to tackle the second order nonlinear 
ODEs provide first integrals only. In this work, we have proposed a novel method to 
identify transformation coordinates which can be used to rewrite the first integral with- 
out explicit time dependence so that it can be integrated by quadrature. Interestingly, 
we showed that the transformation coordinates can be constructed from the first inte- 
gral itself. With this choice, the above procedure can be augmented with other existing 
methods to obtain a complete solution for a given problem. We have also shown the 
applicability of our method in constructing linearizing transformations as well in a sim- 
ple and straightforward way. Finally, we mention that the proposed method can also be 
used to study certain coupled nonlinear oscillators, the details of which will be presented 
separately [Oj. 
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Government of India, sponsored research project. 
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